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ABSTRACT 

The ambient interstellar environment of wind- and supernova-driven superbubbles strongly affects 
their evolution, but its properties are rarely well-determined. We have therefore obtained H I aperture 
synthesis imaging of the environment around three similar, optically-selected superbubble nebulae in 
the Large Magellanic Cloud. The resulting H I maps show that the ambient gas distribution around 
these superbubbles differ to an extreme: DEM L25 shows no neutral shell component, but is nestled 
within an H I hole; DEM L50 shows a massive neutral shell component, but is otherwise within an H I 
void; and DEM L301 shows no correspondence at all between the optical nebula and H I distribution. 
There is also poor correspondence between the H I and optical kinematics. These results strongly caution 
against inferring properties of the ambient neutral environment of individual superbubbles without direct 
observations. Finally, all three objects show some evidence of shock activity. 

Subject headings: ISM: bubbles — H II regions ISM: kinematics and dynamics — ISM: structure 
supernova remnants — Magellanic Clouds 


1. INTRODUCTION 

The evolution of superbubbles created by the super¬ 
novae (SNe) and stellar winds of OB associations is a crit¬ 
ical determinant of the structure and energetics of the in¬ 
terstellar medium (ISM). The coronal component of the 
ISM, or hot ionized medium (HIM) is believed to orig¬ 
inate within such superbubbles and supernova remnants 
(SNRs; e.g., Cox & Smith 1974), hence the evolution of 
these structures determines whether and how this hot gas 
is released into the diffuse HIM. Models for the ISM are 
strongly dependent on whether the HIM is a pervasive, 
dominant component, as envisioned by e.g., McKee & Os- 
triker (1977), or whether it plays a lesser role as favored by, 
e.g., Slavin & Cox (1993), returning to a paradigm closer to 
the two-phase ISM of Field, Goldsmith, & Habing (1969). 
Superbubble evolution also directly impacts the relative 
volumes of enriched, coronal gas in galactic disks vs. ha¬ 
los (e.g., Shapiro & Field 1976). Similarly, superbubbles 
and their related filaments are thought to provide much 
of the structure for the diffuse ISM, including warm and 
cold H I and the warm ionized medium (e.g., Kennicutt 
etal. 1995). Superbubble activity is therefore fundamen¬ 
tal to the ISM phase balance, structure, and kinematics; 
and it is also a principal driver of galactic evolutionary 
processes. 

The standard model for the evolution of these superbub¬ 
bles is an adiabatic, pressure-driven shell with continuous 


mechanical energy injection (e.g., Pikel’ner 1968; Castor, 
McCray, & Weaver 1975). The mechanical power L is pro¬ 
vided by stellar winds and SNe of the parent OB associa¬ 
tion. In this model, an outer shock sweeps up ambient ISM 
into a thin cooled shell, while an inner, reverse shock ther- 
malizes the mechanical energy, thereby generating the hot 
interior whose pressure drives the shell growth. It is gen¬ 
erally assumed that successive SNe can be approximated 
as a continuous energy injection, and that all the available 
SN energy is thermalized (e.g., Mac Low & McCray 1988). 
The standard, adiabatic model then yields simple analytic 
relations governing the shell evolution: 

R oc (L/n) 1/5 f 3/5 , (1) 

v = R<x (L/n) 1/5 t~ 2/5 , (2) 

where R and v are the shell radius and expansion veloc¬ 
ity, respectively, n is the ambient density, and t is elapsed 
time. 

There are several lines of evidence supporting this stan¬ 
dard model for superbubble evolution (see Oey 1999): 
a few young, stellar wind-dominated objects have well- 
constrained L from classification of the stellar populations, 
and show reasonable consistency with the predicted kine¬ 
matics (Oey 1996b). In addition, multiwavelength obser¬ 
vations of highly active star-forming regions like 30 Do- 
radus (Wang 1999) and the starburst galaxy NGC 5253 
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(Strickland & Stevens 1999) show hot gas morphologically 
contained within cooler shell structures. Higher ionization 
species like C IV and Si IV have been detected in lines 
of sight through all superbubbles with available observa¬ 
tions (Chu etal. 1995), suggesting interface zones between 
hot and cold gas. In addition, the superbubble size dis¬ 
tribution derived from the standard model is in remark¬ 
able agreement with observed size distributions of H I 
shells in the Magellanic Clouds (Oey & Clarke 1997; Kim 
etal. 1999). 

However, other evidence suggests that superbubble evo¬ 
lution is often not nearly as simple as described by the 
standard model. The first problem is that the model over¬ 
estimates the shell growth rate for young, wind-dominated 
superbubble nebulae (Oey 1996b; Garcfa-Segura & Mac 
Low 1995; Drissen etal. 1995). This could be due to a 
systematic understimate in the ambient densities n by up 
to a factor of ~ 10, or to a similar overestimate in the 
input power L predicted by the assumed stellar mass-loss 
rates. Secondly, Oey (1996b) identified a class of super¬ 
bubbles that exhibit expansion velocities that are several 
times higher than predicted, relative to their radii. Most of 
these also show anomalously bright X-ray emission, which 
suggests shell acceleration by SNR impacts (Chu & Mac 
Low 1990; Oey 1996b). In such cases, not all the avail¬ 
able SN energy will be thermalized to power the shell ex¬ 
pansion, thus the effect of discrete SNe on shell evolution 
remains unclear. 

Clarifying the structure of the ambient medium is 
clearly essential in determining the origin of these dis¬ 
crepancies. Since most estimates of the ambient n for 
young nebular shells are based on the Ha electron den¬ 
sities, the presence of a massive H I shell component could 
imply a significant underestimate in assumed n. Another 
possibility is that the structure of the ambient medium is 
highly irregular, providing openings to a network of low- 
density channels. In that case, the superbubble could be 
ruptured in one or more locations, depressurizing the inte¬ 
rior and releasing hot gas into these channels. Soft X-ray 
imaging with Einstein and ROSAT suggests this process 
in at least one object, DEM L152 (Chu etal. 1993; Mag- 
nier etal. 1995) in the Large Magellanic Cloud (LMC). 
These minor blowouts would then cause accelerated, high- 
velocity kinematics in excess of the standard model predic¬ 
tions (e.g., Oey & Smedley 1998; Mac Low et al. 1998). Ex¬ 
pansion into the density gradient perpendicular to galac¬ 
tic disks offers an additional systematic explanation for 
anomalously fast shell expansion in face-on galaxies (Silich 
& Franco 1999). 

2. OBSERVATIONS 
2.1. Optical narrowband imaging 

In order to examine their interstellar environment in 
more detail, we therefore obtained H I imaging around 
three nebular superbubbles in the LMC. All three objects 
are among those found by Oey (1996b) to show anoma¬ 
lously high expansion velocities relative to their radii. The 
characteristics of these objects are summarized in Table 0 
columns 1 and 2 give their identifications in the Ha cat¬ 
alogs of Davies, Elliott, & Meaburn (1976) and Henize 
(1956), respectively; and column 3 gives the shell radius 
R from Oey (1996b), assuming a 50 kpc distance to the 


LMC. All three objects contain OB associations in pro¬ 
jection, and their stellar content was classified by Oey 
(1996a). The stars with the earliest spectral types within 
each superbubble are listed in column 4. Finally, col¬ 
umn 5 lists the ratio of predicted to observed Ha luminos¬ 
ity Ana (pred) /Ln a (obs), determined by Oey & Kennicutt 
(1997) from the observed stellar population. 

We obtained optical narrowband imaging of these ob¬ 
jects in the light of Ha, [O III]A5007, [S II]AA6717,6731, 
and also red and green continuum. These observations 
were obtained from the Magellanic Clouds Emission-Line 
Survey (MCELS; Smith et al. 1999), which is nearing com¬ 
pletion at the CTIO/U. Michigan 0.9-m Curtis Schmidt 
Telescope. The observations for DEM L301, which are also 
part of this survey, were published by Oey etal. (2000), 
and similar narrowband imaging was reported by Skel¬ 
ton etal. (1999). The filter bandwidths for Ha, [O III], 
and [S II] are 30A, 40A, and 50A, respectively, with wider 
continuum filters (see Smith etal. 1999); the spatial res¬ 
olution is about 3 — 4", with 2.3" px -1 . Figure [I] shows 
color composite images constructed from the continuum- 
subtracted optical narrowband observations. Ha is shown 
in red, [O III] in green, and [S II] in blue. The nebular 
gas in all three objects is well-defined, highly circular, and 
filamentary. The MCELS narrowband frames of DEM L25 
and DEM L50 were flux-calibrated with two longslit spec¬ 
tra of positions in DEM L25. These were obtained at the 
CTIO 1.5-m telescope on 2001 August 1 under photomet¬ 
ric conditions. 

2.2. HI and continuum observations 

We used the Australia Telescope Compact Array 
(ATCA) to construct aperture synthesis imaging of these 
regions in the 21-cm H I line. As seen in Figure |, DEM 
L25 and L50 are separated by only ~ 0.5°, so for these 
targets we mosaicked a single field from 13 pointings. Our 
field for DEM L301 comprises 7 pointings. We obtained 
the observations in four configurations: 750A (1996 Nov 
27), 750D (1997 Jan 5-7), 1.5C (1997 Sep 3), and 1.5D 
(1997 Mar 26-30). Excluding the data from the baselines 
involving antenna 6, these configurations yield 40 indepen¬ 
dent baselines ranging from 30 to 1470 nr. The 21-cnr line 
observations were carried out for an LMC heliocentric ra¬ 
dial velocity r>hei ~ 297 kms -1 . The bandwidth of 4 MHz 
was divided into 1024 channels, which were then binned 
to 512 channels. These correspond to a velocity coverage 
—33 < Uhei < 627 kms -1 , with resolution 1.65 knrs -1 
and sensitivity of 2.2 K. We simultaneously observed the 
20-cm, 1.380 GHz continuum emission with a bandwidth 
of 128 MHz and full polarization information. 

We observed the primary flux calibrator PKS B1934- 
638 at the start and end of each observing run. For phase 
and amplitude calibration, we observed PKS B0454-810 
and PKS B0407-658. To increase the uv coverage in the 
shortest baselines, we also merged our data with fields from 
the LMC ATCA survey of Kim et al. (1998) and the Parkes 
single-dish survey (Staveley-Smith etal. 2001). For DEM 
L25 and L50, we included 50 fields from the ATCA survey, 
and for DEM L301, we included 37 fields. The pointing po¬ 
sitions of our new data were adopted to coincide with the 
central pointings from the survey data. Altogether, these 
observations yield H I imaging at 50" resolution within 
roughly 40' of each object. 
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Table 1 

Superbubble sample 


DEM 

Henize 

-R(pc) a 

Sp. Type b 

L H a(pred)/L Ha (obs) c 

DEM L25 

N185 

43 

09 V 

8.6 

DEM L50 

N186 CE 

50 

06.5 V((f)) 

2.0 

DEM L301 

N70 

53 

03 If* 

0.2 


a From Oey (1996b). 
b From Oey (1996a). 
c From Oey & Kennicutt (1997). 


The calibration and mosaicking of the radio data were 
carried out with Miriad software. Our data were trans¬ 
formed to the image plane using superuniform weighting 
and a Gaussian taper of 40". We deconvolved the beam 
using a maximum entropy method, and this was then 
merged with the ATCA data and regridded Parkes data. 
We adopted a flux ratio between the Parkes data and our 
data of 1.2 for DEM L25 and L50, and 1.0 for DEM L301. 

3. RESULTS 

The optical emission-line properties of DEM L301 are 
discussed in detail by Oey etal. (2000), who also present 
the individual narrowband images of that object used here. 
Similar CCD narrowband imaging of DEM L301 was ob¬ 
tained and discussed by Skelton etal. (1999). In addition, 
Oey et al. (2000) obtained longslit spectrophotometry and 
compared these with tailored photoionization and shock 
models, showing that the observed line ratios appear to 
result from a combination of shock and density-bounded 
photoionization. 

In contrast, DEM L25 and L50 emit substantially more 
Ha luminosity than apparently can be provided by the 
stars (Table pj). This implies an additional ionization 
mechanism, most plausibly shock heating. Figure jin 
shows a dramatic contrast between DEM L25 and L50 
in [O III] emission: while DEM L25 shows a ratio of 
[O III]A5007/Ha ~ 0.8, the main shell of DEM L50 shows 
almost no detection in [O III]. There is some [O III] emis¬ 
sion from the SNR N186 D, which is the small loop at 
the north end of DEM L50 (see §3.2 below), but Lasker 
(1977) cites the lack of [O III] emission in the main shell 
as evidence that N186 D is physically separate from the 
remainder of the system. The [O III] emission from DEM 
L25 is surprising in view of the stellar populations: DEM 
L25 has only a single 09 V star, whereas DEM L50 has 
one 06.5 V star and several late-type O stars. We would 
have expected [O III] emission from both nebulae, since it 
is even seen in H II regions dominated by 09 stars (Hunter 
1992; Hunter & Massey 1990). However, the extremely low 
ionization parameter in these shell structures is likely to 
quench the photoionized [O III] emission in these objects 
(e.g., Oey etal. 2000; Garnett & Kennicutt 1994). The 
[O III] from DEM L25 is therefore most likely caused by 
an additional heating mechanism, probably shock excita¬ 


tion. 

The nebular kinematics of all three superbubbles have 
been studied extensively (e.g., Rosado etal. 1981; Lasker 
1981; Dopita etal. 1981; Rosado etal. 1982; Meaburn, 
McGee, & Newton 1984; Rosado etal. 1990). All three 
show highly complex nebular velocity structure, and do 
not exhibit simple uniform expansion. The complex kine¬ 
matics were cited as evidence consistent with the sugges¬ 
tion of internal SNR impacts to explain the anomalously 
high velocities (Oey 1996b). 

Figure 0 shows the wide-angle maps of H I column den¬ 
sity 7V(H I) around these regions. The circles indicate the 
positions of the superbubbles at their approximate angular 
size. Figure |J gives a closer look with composite images of 
the optical and H I data. Red is Ha, green is [O III], and 
blue is now H I. It is apparent that the neutral environment 
varies dramatically among the three objects. DEM L25 
and L50 appear to be situated in a large H I void, roughly 
1.4 kpc in diameter, which apparently corresponds to the 
supergiant shell LMC 8 (Meaburn 1980). This structure 
is also clearly seen in the full H I map of the LMC con¬ 
tracted by Staveley-Smith etal. (2001), which combines 
data from the Kim etal. (1998) ATCA survey and the 
Parkes multibeam receiver data. Within this void, there 
is a patch of H I emission corresponding to DEM L50, but 
no apparent H I associated with DEM L25. Figure [| shows 
H I maps of the region over six velocity intervals, which 
are the average over four channels. The mean velocity is 
shown in the upper left of each panel. The environment of 
DEM L301 (Figure pp) shows neutral gas throughout the 
region in projection, out our data again show no apparent 
correspondence in the H I distribution with the nebular 
superbubble (see below). Figure || shows H I maps of this 
region averaged over five velocity channels. Hence, despite 
the similar appearance of the three shells in Ha, their H I 
environments are different to an extreme. 

3.1. DEM L25 

As mentioned above, Oey & Kennicutt (1997) found 
that both DEM L25 and DEM L50 show more Ha emis¬ 
sion than can be accounted for by the late-type O stars 
present within these objects. Both nebulae might there¬ 
fore be expected to be radiation-bounded and associated 
with neutral gas. The environment of DEM L25 is con- 
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sistent with the superbubble being nestled within a par¬ 
tial H I hole (Figure ||), within a region where typically 
N (H I)~ 1.0 x 10 20 cm -2 . While some of the emission 
seen to the NE of DEM L25 probably belongs to an unas¬ 
sociated velocity component, the individual channel maps 
do remain consistent with general H I anticorrelation for 
the superbubble (Figure ||). 

The H I scale height h of the LMC is dependent on ra¬ 
dius, and difficult to determine because of the galaxy’s low 
inclination angle. Kim etal. (1998) estimated a character¬ 
istic h ~ 600 pc, which is weighted toward outer radii and 
assumes that the stars and gas are coupled. A more clas¬ 
sical analysis by Kim etal. (1999) suggests h ~ 180 pc, 
assuming stars and gas are not coupled. Elmegreen, Kim, 
& Staveley-Smith (2001) note a break in the H I power 
spectrum around a scale of 100 pc, which may also indi¬ 
cate a low value of h. Here we adopt a value of h = 300 
pc as a compromise estimate. If the value of N (H I) found 
in the DEM L25 environment is characteristic of the ISM 
in this region, our adopted h implies an ambient volume 
density of n(Hi) ~ 0.05 cm -3 . This low value is related 
to the location of DEM L25 within the supergiant shell 
environment mentioned above. 

There is no emission that unambiguously could corre¬ 
spond to a neutral shell of swept-up material, thus it may 
be that the entire shell is ionized. As seen in Figure [y, 
DEM L25 shows the highest excitation among the three 
objects, with [O III]A5007/Ha ~ 0.8, despite the fact that 
DEM L25 has the latest-type OB stars of the three neb¬ 
ulae. The ratio of observed to predicted Ha luminosity 
is so high, LHa(obs)/THa(pi'ed) = 8.6 (Table [lj), that an 
additional ionization mechanism is unavoidable. This is 
consistent with the high [O III]/Ha ratio, and it is there¬ 
fore less surprisingto find no neutral gas component. 

Figures |p and || show that on the west side of DEM 
L25, the superbubble is interacting with an H I cloud. 
This is apparent both morphologically and from the en¬ 
hanced optical nebular emission, indicating a density in¬ 
crease in the optical shell as would be expected from a col¬ 
lision. The [S II]/H a ratio is also higher in the interaction 
zone, ranging around ~ 1.1, as would be expected from a 
ram pressure-induced shock. The mass of the H I cloud is 
~1.3xl0 4 Mg, as measured within a 230" box centered 
at 04 53 49.4, - 70 00 12 (J2000). Figure @ shows an E-W 
position-velocity (PV) slice across the cloud, showing es¬ 
sentially a single component with a systemic velocity of 240 
kms -1 . The velocity dispersion at the interaction zone is 
about twice as large as for the remainder of the cloud, as is 
qualitatively expected from a collision. The surrounding 
H I spatial distribution suggests that the shell is expanding 
into the cloud, or rather than that the entire superbubble 
has a systemic velocity in collision with the cloud. The 
standard, wind-driven evolution predicts a shell expansion 
velocity v s ~ 15 kms -1 , based on the observed stellar 
population and shell size (Oey 1996b). This value rep¬ 
resents a rough lower limit to the collision velocity, since 
Ha kinematic observations by Rosado et al. (1982) show 
much higher velocity structures that may imply v s up to 
70 kms -1 . We note that no kinematic features are seen in 
the H I observations at these high velocities, again arguing 
that the shell itself is fully ionized. 

To the N and NE of DEM L25 is denser H I gas (Fig¬ 
ure which is a southern extension from a major star¬ 


forming region with high-density gas about 300 pc to the 
north (Figures |2]a and |(i|). Figure |] shows the PV dia¬ 
gram for a N-S slice at position 04 54 28.9, -69 56 21.0, 
to the NE of the shell. Two distinct kinematic compo¬ 
nents are clearly apparent around 245 and 270 kms -1 , 
which undoubtedly correspond to those seen in the unre¬ 
solved observations by Meaburn etal. (1984). These two 
components are continuously present over most of the H I 
emission in the large mass immediately to the NE of DEM 
L25 and contrast with the single component of the cloud 
to the W. The existence of two, continuous, distinct veloc¬ 
ity components is puzzling, but is a common feature that 
is also seen, for example, in the Small Magellanic Cloud 
(Staveley-Smith etal. 1997). 

3.2. DEM L50 

DEM L50 shows an H I configuration that is virtually 
inverse to that of DEM L25. DEM L50 clearly shows as¬ 
sociated H I, plausibly in a neutral shell envelope outside 
the nebula, while the surrounding environment has a low 
density. For a single 1.65 kms -1 channel, our 3er limit 
is 2.0xl0 19 cm -2 , although the effective detection limit 
depends on the kinematic range of the structure. Binning 
over four channels (Av = 6.6 kms -1 ) as shown in, e.g., 
Figured gives a nominal detection limit of4.0xl0 19 cm -2 . 
This limiting ambient H I column density corresponds to 
a volume density of n(Hl) <, 0.02 cm -3 , an extremely 
low value that is again associated with the interior of the 
supergiant shell environment. Figure d shows the excel¬ 
lent spatial correspondence of the Ha (image) and N (H I) 
(contours) maps. The total integrated emission of 282 Jy 
kms -1 within a 610" box implies a total H I mass of 
1.7x10 5 Mq for this material. For comparison, the ambi¬ 
ent density of 1.4 cm -3 estimated by Oey (1996b) and 50 
pc radius imply a swept-up mass of 1.6 xlO 4 M 0 . This 
therefore suggests that the ambient density was underes¬ 
timated by an order of magnitude from the Ha emission. 
Alternatively, it may that most of the H I originates in a 
structure exterior to the superbubble, although this seems 
unlikely because of the spatial correlation between the Ha 
and H I emission in Figure |j. 

Figure [To] shows PV diagrams for E-W and N-S slices 
across the center of DEM L50. Most of the material has 
a velocity ~ 247 kms -1 . The N-S slice (Figure [To|a,) 
shows a gradient to ~ 260 kms -1 near the edges, sug¬ 
gesting that most of the central shell material may cor¬ 
respond to a near-side shell expansion velocity of v s ~ 
13 kms -1 . Based on the observed stellar population 
and shell radius, the standard wind-driven evolution pre¬ 
dicts v s ~ 20 kms -1 for the originally assumed ambient 
n = 1.4 cm -3 , and v s ~ 16 kms -1 for n increased by 
an order of magnitude (Oey 1996b). Both values are con¬ 
sistent with a near-side H I expansion. The Ha veloc¬ 
ity structure is again complex (Rosado etal. 1990); while 
the reported values are broadly consistent with the prin¬ 
cipal H I component around ~ 240 kms -1 , it is difficult 
to identify the optical kinematic structures with the H I. 
Rosado et al. suggest a possible optical shell expansion of 
v 3 ~ 25 ± 9 kms -1 , reasonably consistent with the ten¬ 
tative near-side H I expansion. Howver, it is also possible 
that the entire superbubble may have a systemic velocity 
offset from the ambient ISM, or that the shell kinematics 
are more complicated. 
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Figure |lT] shows a PV diagram for the SNR N186 D, 
obtained from a N-S slice at position 05 00 12.3, 70 06 
22. Here we do clearly see a ~ 20 kms -1 line split¬ 
ting over a region extending about 4'. However, this is 
puzzlingly inconsistent with Ha observations suggesting 
v s ~ 90 kms -1 for the SNR, found in the Ha Fabry-Perot 
data of Laval et al. (1989). Their data generally appear to 
show two components around 160 and 320 kms -1 , and 
there are no reported velocity structures in the range 240 
to 270 kms -1 , which are seen in our H I observations. Our 
data cube has limiting extremes of 190 and 336 kms -1 , 
so we cannot adequately probe the same kinematic range 
as the optical data, but we detect no structures near these 
extremes in our cube, within our effective detection limit 
of ~ 4 x 10 19 cm -2 . 

3.3. DEML301 

The H I environment of DEM L301 is still different in 
character from either that of DEM L25 or DEM L50. As 
seen in Figure [l^, there is no apparent correspondence at 
all between the H I emission and nebular emission, whether 
as a correlation or anticorrelation. Nor can we identify 
any correspondence in the individual channel maps, (Fig¬ 
ure ||). Figure |l3| shows PV diagrams for E-W and N-S 
slices across the center of the optical nebula. There is 
no kinematic structure readily identifiable with the super¬ 
bubble either. Optically, DEM L301 again shows complex, 
high-velocity components in Ha, with velocity differences 
up to 100 kms -1 (Rosado etal. 1981). Since the object is 
density-bounded and leaking ionizing photons (Table Q), 
the shell is probably fully ionized, which helps explain the 
lack of corresponding H I structure. However, this would 
still imply that the H I distribution is clumpy or patchy, 
to explain the absence of an H I hole. Typically in this re¬ 
gion N(H. I)~ 1.3 x 10 21 cm -2 , suggesting a mean volume 
density of 0.7 cm -3 . 

4. 20-CM CONTINUUM OBSERVATIONS 

Our 20-crn (1380 MHz) continuum observations of the 
three superbubbles show excellent spatial correspondence 
with the Ha distributions, as expected for the nebular 
thermal continuum. We measure flux densities within our 
128 MHz bandwidth of 84 ± 17 mJy and 173 ± 35 mJy for 
DEM L25 and DEM L301, respectively, excluding point 
sources. The flux density of DEM L50, excluding one point 
source, is 229 ± 46 mJy, of which 99 ± 20 mJy corresponds 
to the SNR in the northern region. 

Since all three superbubbles have been suggested to be 
affected by shocks, there may be a non-negligible contribu¬ 
tion to the continuum from non-thermal emission. To eval¬ 
uate this possibility, we constructed spectral index maps 
using our 1380 MHz data and 843 MHz data from the Mo- 
longlo Observatory Synthesis Telescope (MOST) survey of 
the Magellanic Clouds (Turtle etal. 1998). The results are 
noisy, but broadly suggest that DEM L25 and the SNR 
N186 D in DEM L50 have spectral indices roughly around 
a ~ —0.8 for flux density S oc v a : whereas DEM L301 
and the remainder of DEM L50 have lower spectral in¬ 
dices a ~ —0.3. This is consistent with the stronger shock 
signatures seen for DEM L25 and SNR N186 D, specifi¬ 
cally, the higher [O III]/Ha, [S II]/Ha, overproduction of 
Ha, and higher-velocity kinematics. The data thus quali¬ 
tatively suggest the presence of some non-thermal emission 


for at least DEM L25 and SNR N186 D. These results are 
consistent with the spectral indices reported by Filipovic 
etal. (1998) in their multi-band radio continuum survey. 
For DEM L25 and DEM L50 (unresolved), they obtain 
a = —0.67 ± 0.09 and a = —0.73 ± 0.19, respectively, 
while for DEM L301 they find a = —0.02 ± 0.03. These 
values may be compared to the mean a = —0.15, standard 
deviation a = 0.31 for H II regions; and mean a = —0.43, 
a = 0.19 for SNRs, reported by the Filipovic etal. (1998) 
survey. 

5. CONCLUSION 

We have obtained H I observations of the environ¬ 
ment around three optically-selected superbubbles in the 
LMC. The optical nebulae have similar properties: all are 
roughly 50 pc in radius and show optical kinematic struc¬ 
tures with anomalously high expansion velocities that con¬ 
trast with values expected from the standard evolutionary 
model for wind and SN-generated superbubbles. We find 
that the H I environment within roughly 500 pc of the 
objects vary to an extreme. DEM L25 shows no H I asso¬ 
ciated with the superbubble itself, but appears to be nes¬ 
tled within an H I hole, and the shell is interacting with 
an adjacent H I cloud; DEM L50 shows an H I component 
to the shell that is ~ 10 times more massive than the ion¬ 
ized shell, but is otherwise in a large H I void; and DEM 
L301 shows no correlation of any kind between the opti¬ 
cal and H I distributions. Furthermore, the high-velocity 
kinematics observed optically for all three objects are not 
seen in any of our H I data. These findings strongly cau¬ 
tion against inferring properties of the neutral interstellar 
environment around any individual superbubbles without 
direct observations. 

All three superbubbles have reported signatures sugges¬ 
tive of shock activity, including supersonic nebular veloci¬ 
ties. DEM L25 and DEM L50 show elevated Ha emission 
in excess of that predicted from photoionization by their 
classified stellar populations, and DEM L50 encompasses 
the known SNR N186 D on its northern edge. DEM L25 
also shows elevated [O III] /Ha and all objects show high 
[S II]/Ha. The emission-line spectrum of DEM L301 was 
modeled in detail by Oey etal. (2001), who found evi¬ 
dence of contributions by shock excitation. We examined 
our 1380 MHz (20 cm) continuum in conjunction with 843 
MHz continuum observations from the MOST survey to 
obtain rough, spatially-resolved estimates of the spectral 
indices. These suggest some presence of a non-thermal 
component for at least DEM L25 and the SNR N186 D. 

We mentioned above that the new H I data do not 
show the complex, high-velocity shell components that 
are evident in the optical nebular data for these objects. 
However, the uneven H I distributions for DEM L25 and 
DEM L301 are consistent with the high ionization state 
of these nebulae, suggesting that these shells are fully ion¬ 
ized. Their patchy H I remains consistent with the optical 
filaments corresponding to incipient blowout features or 
material accelerated by SNR impacts, as has been previ¬ 
ously suggested. The lack of correspondence between the 
nebular and H I kinematics is more problematic for DEM 
L50, although the dominant Ha and H I components are 
broadly consistent with each other. The H I velocity pro¬ 
file suggests that the neutral component may correspond 
to the approaching side of the shell. 
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Fig. 1.— Composite narrow-band images of a) DEM L25 and DEM L50; and 6) DEM L301. Ha is shown in green, [O III] in blue, and [S II] 
in red. Figure LL is constructed from the continuum-subtracted images with arbitrary intensity scaling. Beware that the relative intensities of 
the three bands vary between the two fields. North is up, east to the left; see Figures pi pi y, and fi2lfor angular scale. 


Fig. 2.— H I column density maps of the environment around a) DEM L25 and DEM L50; and b ) DEM L301. The circles roughly indicate 
the position and size of the objects, with DEM L25 on the west side of Figure Ha. 


Fig. 3.— Composite images of a) DEM L25 and b) DEM L50, with Ha shown in red, [O III] in green, and H I in blue, using arbitrary 
intensity scaling. 


Fig. 4.— H I maps of the region around DEM L25 and L50. Each panel shows the average of four channels, whose mean heliocentric 
velocity in kms —1 is shown in the upper left. The circles indicate the positions of DEM L25 and L50, with the former on the west side. 


Fig. 5.— H I maps for the region around DEM L301. Each panel shows the average of five velocity channels. 


Fig. 6.— 7V(H I) contour map of DEM L25, superimposed on Ha. Contour intervals are 2 x 10 20 cm 2 , beginning at 2 x 10 19 cm 2 . 


Fig. 7.— Position-velocity diagram oriented E-W across the cloud to the W of DEM L25. The center is at 04 53 14.7, -69 59 35, with E 
to the left (negative offsets). 


Fig. 8.— Position-velocity diagram for the H I mass NE of DEM L25, centered at 04 54 29.9, -69 55 35 (J2000). N is to the left (negative 
offsets). 


Fig. 9.— 7V(H I) contour map of DEM L50, superimposed on Ha. Contours are as in Figure ( 


Fig. 10.— Position-velocity diagrams across the center of DEM L50, oriented a) E-W; and b) N-S. E and N are to the left (negative offsets). 


Fig. 11.— Position-velocity diagram across the SNR N186 D. The slice is oriented N-S, centered at 05 00 12.3, -70 06 23 (J2000), with N 
shown to the left. 


Fig. 12.— N(H I) contour map of DEM L301, superimposed on Ha. Contour intervals are 3 x 10 2 ° cm 2 , beginning at 2 x 10 19 cm 2 . 


Fig. 13.— Position-velocity diagram across the center of DEM L301, oriented a) E-W; and b) N-S. E and N are shown to the left. 
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